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Abstract: We perform a comprehensive study of the impact of new-physics opera- 
tors with different Lorentz structures on decays involving the b — > sfi + [i~ transition. 
We examine the effects of new vector-axial vector (VA), scalar-pseudoscalar (SP) "' 
and tensor (T) interactions on the differential branching ratios and forward-backward 
asymmetries (A F b's) of B° s ->■ B Q d ->■ X s /i + /i", B° s ->■ /i + /U~7, B° d ->■ K^pT , 

and Bj — > K*fi + fi~ , taking the new-physics couplings to be real. In B® — > K* n + n~ , 

(2) 

we further explore the polarization fraction Jl, the angular asymmetry A T , and 
the longitudinal-transverse asymmetry Alt- We identify the Lorentz structures that 
would significantly impact these observables, providing analytical arguments in terms 
of the contributions from the individual operators and their interference terms. In 
particular, we show that while the new VA operators can significantly enhance most 
of the asymmetries beyond the Standard Model predictions, the SP and T operators 
can do this only for Afb in B® — > K[i + [i~ . 
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1. Introduction 



In recent years, there have been quite a few measurements of quantities in B decays 
which differ from the predictions of the Standard Model (SM) by ~ 2a. For example, 
in B — > 7iK, the SM has some difficulty in accounting for all the experimental 
measurements [|I|]. The measured indirect (mixing- induced) CP asymmetry in some 
b — > s penguin decays is found not to be identical to that in B® — > J /tpK s ||, [3|, f|, 
counter to the expectations of the SM. While the SM predicts that the indirect CP 
asymmetry in B® — > J/ip4> should be ~ 0, the measurement of this quantity by 
the CDF and D0 collaborations shows a deviation from the SM ||. One naively 
expects the ratio of transverse and longitudinal polarizations of the decay products 
in B -> (f>K* to be f T /f L < 1, but it is observed that f T /f L ~ 1 |, |. It may 
be possible to explain this value of f T / II within the SM, but this is not certain. 
Finally, the recent observation of the anomalous dimuon charge asymmetry by the 
D0 collaboration M also points towards some new physics in B s mixing that affects 
the lifetime difference and mixing phase involved therein (for example, see Ref. 
Though none of the measurements above show a strong enough deviation from the 
SM to claim positive evidence for new physics (NP), they are intriguing since (i) the 
effects are seen in several different B decay channels, (ii) use a number of independent 
observables, and (iii) all involve h — > s transitions. 

A further hint has recently been seen in the leptonic decay channel: in the 
exclusive decay B® — > K* n + ji~ , the forward-backward asymmetry (A FB ) has been 
found to deviate somewhat from the predictions of the SM O, [K 



|1J. This 

is interesting since it is a CP-conserving process, whereas most of the other effects 
involve CP violation. Motivated by this tantalizing hint of NP in B® — > K*{i + fi~ , we 
explore the consequences of such NP in related decays. We do not restrict ourselves to 
any particular model, but work in the framework of effective operators with different 
Lorentz structures. 

If NP affects B® — > K*fi + fi~, it must be present in the decay b — > s/i + /i~, and 



will affect the related decays B® — > , B® —> X s fi + fi , B® — > ji 



H 7. 



and 59 



Kn + n . The analyses of these decays in the context of the SM as well as in some NP 
models have been performed in the literature: B° s — > fj, + /j, 



15, m, 17, 18, 19, 20, 21 



22], H, H H, B° d -> Xsii+iT ]26 



9} gg, 0, ||§§ |35|, B° s fi+fi'j 

0, n, gg, 15 gg, n, m, 5° Rfi+fi- n, n, n, 0, m, 0, 5^ 

K*fi+fi- pa, m, m m, et\, m m, pm m 



I B3L B7|. Correlations 

between some of these modes have been studied in Refs. p?| ^DL [7D| . 

In this paper, we consider the addition of NP vector-axial vector (VA), scalar- 
pseudoscalar (SP), and tensor (T) operators that contribute to b — > sfi + fi~, and 
compute their effects on the above decays. Our aim here is not to obtain precise 
predictions, but rather to obtain an understanding of how the NP affects the observ- 
ables, and to establish which Lorentz structure(s) can provide large deviations from 
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the SM predictions. Some of these effects have already been examined by some of us: 
for example, new VA and SP operators in B° s — > ji + n~ new VA and SP operators 
in — > /i + /i~7 the correlation between B® — > and 5° — > Kfi + fi~ with 

SP operators 7C], large forward-backward asymmetry in B® — >■ Kl^lC from 
T operators P9[ , and the contribution of all Lorentz structures to 5° — > K*n + n~ , 



with a possible explanation of the Aps anomaly p2 |. Here we perform a combined 
study of all of these decay modes with all the Lorentz structures, consolidating and 
updating some of the earlier conclusions, and adding many new results and insights. 
Such a combined analysis, performed here for the first time, is crucial for obtaining a 
consistent picture of the bounds on NP and the possible effect of NP on the observ- 
ables of interest. While observables like the differential branching ratio (DBR) and 
ApB(q 2 ) by themselves are sensitive to NP, we also examine the correlations between 
them in the context of NP Lorentz structures. 

A full angular distribution of B^ — > K* \i~ allows us access to many inde- 
pendent observables, and hence to multiple avenues for probing NP. We present 
here for the first time the full angular distribution, including all the NP Lorentz 
structures, for this decay mode. This leads to the identification of observables that 
could be significantly influenced by specific Lorentz structures of NP. In addition to 
the DBR and Aps-, we also examine the longitudinal polarization fraction fi and 

( 2) n 

the angular asymmetry A T , introduced recently in Ref. pT7| . We further analyze 
the longitudinal-transverse asymmetry Alt, which, as we will argue, has very small 
hadronic uncertainties. 

Hadronic uncertainties often are the main source of error in the calculation of SM 
predictions of a quantity, and make the positive identification of NP rather difficult. 
In this paper, for B® — > K[i + n~ we use the form factors from light-cone sum rules. 
For B® — > K*fi + fi~, we use the form factors obtained from QCD factorization at low 
q 2 , and those from light-cone sum rules at high q 2 . The latest next-to- leading order 
(NLO QCD) corrections |7l] have not been included. These corrections would affect 
the central values of the SM predictions to a small extent, while also decreasing the 
renormalization-scale uncertainty. However, since our primary interest is looking for 
observables for which the NP effects are large, a LO analysis is sufficient at this stage. 
In our figures, we display bands for the SM predictions that include the form-factor 
uncertainties as claimed by the respective authors. 

In addition to the form-factor uncertainties, the SM prediction bands also include 
the uncertainties due to quark masses, Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements and meson decay constants. In our figures, these bands are overlaid with 
some examples of the allowed values of these observables when NP contributions are 
included. This allows the scaling of these uncertainties to be easily visualized. It 
turns out that in many cases, the results with the NP can be significantly different 
from those without the NP, even taking into account inflated values for the hadronic 
uncertainties. We identify and emphasize such observables. We also show that the 
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hadronic uncertainties in several of these observables are under control, especially 
when the invariant mass of the muon pair is small and one can use the limit of large- 
energy effective theory (LEET) . This makes such observables excellent probes of new 
physics. Also, since all the observables are shown as functions of q 2 , we have the 
information not just about the magnitudes of the observables, but also about their 
shape as a function of q 2 , where some of the uncertainties are expected to cancel out. 

In this paper, we restrict ourselves to real values for all the NP couplings, and 
study only the CP-conserving observables 1 . In section Q we examine the various SM 
and NP b — > s/i + /i~ operators, and give the current constraints on the NP couplings. 
The effects of the NP operators on the observables of the decays are discussed in 
the following sections: 5° -)■ jj^yT (Sec. |), B Q d ->• X s fi + fi~ (Sec. g), B° s ->■ ,u + /i~7 
(Sec. ^|), B® — > Kji + n~ (Sec. U), an d — > K* n + ji~ (Sec. |7|). Our notation in these 
sections clearly distinguishes the contributions from VA, SP and T operators and 
their interference terms, which offers many insights into their impact on modifying 
the observables. We give the details of the calculations involved in sections ^|-[7| in 
the appendices 0-0, respectively, for the sake of completeness and in order to have 
a clear consistent notation for this combined analysis. In Sec. |^, we summarize our 
findings and discuss their implications. In particular, we point out the measurements 
which will allow one to distinguish among the different classes of NP operators, and 
thus clearly identify which type of new physics is present. 



2. b —> sfi + fi Operators 

2.1 Standard Model and New Physics: effective Hamiltonians 

Within the SM, the effective Hamiltonian for the quark-level transition b — > sy + yT 
is 

n s eS M = ~£r V*V tb fec^O^) + C 7 -^- 2 [sa, u {m s P L + m b P R )b] 
* i=i 

+ C 9 ^ {*fP L b) w + Cio ^ (sYP L b) } , (2.1) 

where Plr = (1 =F 7s)/2. The operators Oi {i = 1,..6) correspond to the Pi in 



Ref. [31], and = m^/i) is the running 6-quark mass in the MS scheme. We use 
the SM Wilson coefficients as given in Ref. pTf . In the magnetic dipole operator 
with the coefficient CV, we neglect the term proportional to m s . 



1 The CP- violating observables, with complex values of the couplings, are treated in the com- 
panion paper j72[ . 
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The operators Oj, i = 1-6, can contribute indirectly to b — > s{i + fi and their 
effects can be included in an effective Wilson coeficient as [^l[ 



Cf = C 9 (m b ) + h(z, m c ) ( to x + C 2 + 6 C 3 + 60 C 5 

4 64 \ 

^h(z, m b ) [ 7C 3 + ~C 4 + 76C 5 + -CA (2.2) 

-/i(z, 0) C 3 + -C 4 + 16 C 5 + — C 6 + -C 3 + — C 5 + — C 6 . 



9 27 

Here 2; = g 2 /m 2 , and m g = m q /m b for all quarks g. The function h(z,m) represents 
the one-loop correction to the four-quark operators 0\-0§ and is given by [^7, 

h(z, rh) = — In — mm + — + —x (2-3) 

9 9 



~(2 + x)\l -x\ 1/2 




4 




-x 




9 






— iix \ 







for a; < 1 



9 2arctan t= , for a; > 1 , 

1 -'a;—!' 



where x = Am 2 /z. In the numerical analysis, the renormalization scale /i& is varied 
between m b /2 and 2mft. Note that in the high-g 2 region one can perform an operator 
product expansion (OPE) in 1/Q with Q = (rribyq 2 ') [73, fflj . Numerically the 
results of Refs. []73|, [74| differ little from those in Eq. ( |2.2|) and so we use the above 
expression for the entire range of q 2 . An analysis of b — > sfi + fi~ where the OPE in 
the high-g 2 region is used can be found in Refs. |>|, |6[] . 

We now add new physics to the effective Hamiltonian for b — > s/i + /i~, so that it 
becomes 

n cS (b SM + /0 = n'J 1 + h v cS a + n s J + , (2.4) 

where is given by Eq. fl2.1|) , while 

nil = ^ V*V tb {R v {arfP L b) W + &a pry,M 

+ R' v (sYP R b) Jty^n + R' A {s^P R b) /i 7M 75/i} , (2.5) 

U S J = ^ V*V tb {r s (sP R b) ^ + R P (sP R b) wi* 

+ R' s (sP L b) ^ + R'p (sP L b) /i7 5 ^} , (2.6) 
Hj s = ~j!r ^ V*V tb [Crisa^fla^fi + iC TE (sa, u b)rua a ^ e^} (2.7) 

are the new contributions. Here, Rv, Ra, R'v, R'a, R-s, Rp, R'si R'p,Ct and Cte are 
the NP effective couplings. We do not consider NP in the form of the O7 = 
sa a P P R b F a p operator or its chirally-flipped counterpart 0' 7 = sa a/3 Pib F a p. This 
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is because there has been no hint of NP in the radiative decays B — > X s 7, K^-y 
which imposes strong constraints on [Cf |. This by itself does not rule out the 



possibility of a flipped-sign CS S scenario. However this solution can be ruled out at 



3<j from the decay rate of B — > X s £ + £~ if there are no NP effects in Cg and C\q [75 
Thus, NP effects exclusively in C7 cannot provide large deviations from the SM. The 
impact of 0' 7 on the forward-backward asymmetry in B d — > K*fi + ^, together with 
other observables, was studied in Ref. J5CB - 

Note that the operators with coefficients Ry and Ra have the same Lorentz 
structure as those in the SM involving Cg and C w , respectively [see Eq. (p.l|)1, so that 
any measurement will be sensitive only to the combinations (Cg + Ry) or (Cio + Ra)- 
For simplicity, in our numerical analysis of the observables of various decays, these 
couplings are taken to be real. As a consequence, the results in this paper would be 
the same if the corresponding CP-conjugate decays were considered. However, for 
completeness, the expressions allow for a complex-coupling analysis. 

When calculating the transition amplitudes, for the leptonic part we use the 
notation 

L" = (fj,+(p + )ii-(p^\fl^fj,\0) , = ( / i+(p + ) / i-(p„)|^7 5 /i|0) , 

L ee (fi + (p + )^(p^\M0) , L 5 = (//+( P+ K(p_)|/2 T ViO> , (2.8) 

L^ = (//+(p + )/x-(p_)|^>|0) . 

2.2 Constraints on NP couplings 

The constraints on the NP couplings in b — > sfi + fi~ come mainly from the upper 
bound on the branching ratio B(B® — > /U + /i~) and the measurements of the total 
branching ratios B(B° d ->• X s /i+/i~) and B{B° d ->• /?/i + /i") |7|, |77|, |78|, ^ || : 



B(B° S p + p-) < 3.60 x 10" 8 (90% C.L.) , (2.9) 

v d \ (0.44 ±0.12) x 10~ 6 (high g 2 ) ' v ; 

B(B° d ^V) = (4.5±1;=) x HT 7 , (2.11) 

where the low-g 2 and high-g 2 regions correspond to 1 GeV 2 < q 2 < 6 GeV 2 and 
q 2 > 14.4 GeV 2 , respectively, where q 2 is the invariant mass squared of the two 
muons. The constraints from the first two quantities above have been derived in 
Ref. [52] . Here we also include the additional constraints from B(B d — > K^i + ^r). 
The three decays above provide complementary information about the NP operators. 
For the SM predictions here, we use the latest NNLO calculations. Note that the 
measurements for B(B d — > K* / u + ^) are also available [11, 1~2[ . However, the form- 



factor uncertainties in B d — > K*fi + fi~ are rather large, and as a result the constraints 
due to this decay mode are subsumed in those from the other three modes. 

The constraints on the new VA couplings come mainly from B(B d — > X s fi + fi~) 
and B(B d — > K^l + ^t). Their precise values depend on which NP operators are 
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Figure 1: The constraints on the couplings Rv,Ra (left panel) and R' V ,R A (right panel) 
when only primed or unprimed couplings are present. 

assumed to be present. For example, if only Ry,A or only R' VA couplings are present, 
the constraints on these couplings take the form shown in Fig. [I]. For Ry,A, the 
allowed parameter space is the region between two ellipses: 

|itV + 3.6| 2 1^-4.0| 2 \R V + 2.8\ 2 |i?A-4.1| 2 
^ (4.7) 2 (4.8) 2 ' (6.5) 2 (6.6)2 & > ^ ■ J 

while for -Ry^, the allowed region is the intersection of an annulus and a circle: 

22.2 £ + 3.6| 2 + \R' A -4.0| 2 £ 56.6 , l^ 2 + |^| 2 £ 17 . (2.13) 

If both Rv,a and R' V a are present, the constraints on them get individually weakened 
to 

|iV + 2.8| 2 ]^-4.1| 2 

(6.5) 2 + (6.6) 2 ^ ' 1 J 

and 

\R' v \ 2 +\R' A f&40, (2.15) 

respectively 2 . 

For the SP operators, the present upper bound on B(B° S — > /U + /i~) provides the 
limit 

\R S - R' s \ 2 + \Rp ~ R'p? & 0.44 , (2.16) 
where we have used / Bs = (238.8 ±9.5) MeV |I| and IV^V^I = 0.0407± 0.0010 



This constitutes a severe constraint on the NP couplings if only -R^p or R' s P are 
present. However, if both types of operators are present, these bounds can be evaded 



2 Note: the constraints on Rv,a obtained here are milder than those obtained in Ref. |Q using 
B(B® — > (K , K*) /i + pr). This is because Ref. Ji3| had neglected the interference terms between 
the SM and new physics VA operators. Their inclusion relaxes the stringent constraints therein. 
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due to cancellations between the Rs,p and R' s P . In that case, B(B% — >■ X s fi + fi~) and 
B(B® — )■ Kfi + fi~) can still bound these couplings. The stronger bound is obtained 
from the measurement of the latter quantity, which yields 

\Rs\ 2 +\Rp\ 2 &9 , Rs^R's, Rp^R' p . (2.17) 

Finally, the constraints on the NP tensor operators come entirely from B(B® — > 
X s n + When only the T operators are present, 

\C T \ 2 + 4\C TE \ 2 £ 1.0 . (2.18) 



Although the bounds presented in this section for VA, SP and T couplings are 
obtained by taking one kind of Lorentz structure at a time, in our numerical analysis 
for scenarious where we consider combinations of two or more kinds of Lorentz struc- 
tres, we use the allowed parameter space obtained by considering the corresponding 
combined Lorentz structures. 

We now analyze the b — > sfi + fi~ modes in detail and present our results. As 
explained in the Introduction, the figures have the SM prediction bands overlaid 
with the predictions for specific allowed values of NP couplings. The SM band 
is generated by varying the form factors within their ranges as predicted by the 
respective authors, while the CKM matrix elements, quark masses and meson decay 
constants are varied within their 1.6<r allowed values. 



3. B° s -+ fi+fi- 

In this section we examine the NP contributions to B® — > fi + fi~. Within the SM, 
B® — > / u + is chirally suppressed. The SM prediction for the branching ratio is 
B{B° a -> fi + fi-) = (3.35 ± 0.32) x 1(T 9 [f2|. The Tevatron gives an upper bound 
on its branching ratio (BR) of 3.6 x 1(T 8 at 90% C.L. © 0> HI- This deca y can 



be observed at the Tevatron only if NP enhances its BR above 10~ 8 . LHCb is the 
only experiment which will probe B(B® — > // + ^r) down to its SM value. It has the 
potential for a 3a observation (5a discovery) of B® — > with ~ 2ftr x (~ 6fbT 1 ) 

of data [0. LHCb therefore has the potential to observe either an enhancement or 



a suppression of B(B® — > ). It can observe B° s — » as long as its BR is 

above 1.0 x 10~ 9 . 
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3.1 Branching ratio 

The transition amplitude for B° s — > is given by 



-i) 



AG F a err , 

V2 4tf 



(0 \s llll5 b\ B° s (p)) (-Cf -Ra + R' a )L^ 
+ <0 \s l5 b\ B° s (p)) [(Rs - R's)L + (Rp ~ R' P )L 5 } \ , (3.1) 



where L 5fl , L and L 5 are defined in Eq. ( ]2.8|) . Using the matrix elements [15 



<0 \s W b\ B° s (p)) = ip„ f Bs , <0 \s l5 b\ B° s (p)) = -i I 



mb + m s 



the calculation of the BR gives 



Am, 



in 



64vr 3 



Rs — R Q 



\v tb v; s \ 2 



x 



rn 



B s 



Rp — Rn 



2m 



m b + m s 



H t~{Cio + Ra — R'a, 



rn 



(3.2) 



(3.3) 



Clearly, NP in the form of tensor operators does not contribute to B Q S — > fi + fi~. From 
Eq. ( p.3|) and the constraints on NP couplings obtained in Sec. |2.2j , one can study 
the effect of new VA and SP couplings. 

Since the NP contribution from VA operators is suppressed by a factor of ~ 
m^/mb compared to that from the SP operators, the effect of SP operators dominates. 
Both VA and SP operators can suppress B(B° S — > significantly below the SM 

prediction. However while VA operators can only marginally enhance B(B® — > fi + fj>~) 
above 10~ 8 , making the decay accessible at the Tevatron in an optimistic scenario, the 
SP operators can enhance the branching ratio even up to the present experimental 
bound. Indeed, the strongest limit on the SP couplings comes from this decay. This 
strong limit prevents the SP operators from expressing themselves in many other 
observables, as we shall see later in this paper. 

3.2 Muon polarization asymmetry 

The longitudinal polarization asymmetry of muons in B® — > / u + pT is defined as 

N R - N L 



A 



LP 



N R + N L ' 



(3.4) 
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where Nr (Nl) is the number of yU~'s emerging with positive (negative) helicity. Alp 
is a clean observable that is not suppressed by m M /mB a only if the NP contribution 
is in the form of SP operators, such as in an extended Higgs sector. 



A L p for the most general NP is [7C 



A 



LP 



2 L _4m^ Re 
Y m B s 


(Rs-R' s \ 


( Rp : R,p + ^(c w + r a r'S) 

\ mb+rris m B a v A> ) 


(l _ 4m l \ 
\ m B 3 ) 


Rs-R's 

m b +m s 


2 

+ 


R P -R' P + 2m a ( C W + R A R>) 
mt,+m a m B a v A ' 


2 



(3.5) 



From the above equation, we see that Alp can be nonzero if and only if R$ — R'g ^ 
0, i.e. there must be a contribution from NP SP operators. (Within the SM, SP 
couplings are negligibly small, so that A L p — 0.) 

The present upper bound on B(B° S — >■ /i + yU~) puts no constraint on App, and it 



can be as large as 100% [f70" |. A^p can be maximal even if B(B® — >■ ) is close 



to its SM prediction. Therefore, in principle ^4pp can serve as an important tool to 
probe NP of the SP form. However, in order to measure its polarization, the muon 
must decay within the detector. This is not possible due to the long muon lifetime 
(cr for the muon is 659 m). Hence in practice, this quantity is not measurable at 
current detectors. 



4. B° d X s ^ + fi- 

The BR of B® — > X s fi + [i~ in the low-g 2 and high-g 2 regions has been measured to 
be 0, ^ 

'1.49 ± 0.50^) x 10- 6 , (Belle) , 
B(B XJ + £-) low q2 = <( (1.8 ± 0.7 ± 0.5) x 10~ 6 , (BaBar) , (4.1) 

(1.60 ± 0.50) x 10~ 6 , (Average) . 

(0.42±0.12±°$) x 10~ 6 , (Belle) , 
B(B X s £ + £-) high q 2 = <{ (0.50 ± 0.25±°$) x 10~ 6 , (BaBar) , (4.2) 

(0.44 ± 0.12) x 10" 6 , (Average) . 

The SM predictions for B(B — >• X s n + ^r) in the low-g 2 and high-g 2 regions are 



;i.59 ± 0.11) x 10~ 6 and (0.24 ± 0.07) x 10~ 6 , respectively |3 . 

Apart from the measurement of the total BR of B® — > X s n + ^, which has 



already been used to restrict the VA and T operators in Sec. |27|, the differential 
branching ratio (DBR) as a function of q 2 also contains valuable information that 
can help us detect NP. In particular, the SM predicts a positive zero crossing for App 
in B® — y X s fi + fi~ in the low-g 2 region, i.e. for q 2 less than (greater than) the crossing 
point, the value of A FB is negative (positive). This zero crossing is sufficiently away 
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from the charm resonances so that its value can be determined perturbatively to an 
accuracy of ~ 5%. The NNLO prediction ]33| for the zero of App(q 2 ) is (taking 
m b = 4.8 GeV) 

( g 2 ) = (3.5 ±0.12) GeV 2 . (4.3) 

This quantity has not yet been measured. However, estimates show that a precision 
of about 5% could be obtained at a Super- factory p5| . A deviation from the zero 
crossing point predicted above will be a clear signal of NP. 



4.1 Differential branching ratio and forward-backward asymmetry 

After including all the NP interactions, and neglecting terms suppressed by m^/rrib 
and m s /mb, the total differential branching ratio dB/dz is given by 




Total 




B 



SM-VA 



+ By a + Bqp + B- 



(4.4) 



where the quantities B depend on the SM and NP couplings and kinematic vari- 
ables. The complete expressions for these quantities are given in Appendix [A]. The 
subscripts denote the Lorentz structure (s) contributing to that term. 
The forward-backward asymmetry in B® — > X s fi + fi~ is 



A FB {q 2 



Jq d cos 9 



d 2 B 



f 1 dq 2 d cos 9^ 



J_ x d cos 9 



d 2 B 



I 1 dq 2 d cos 9^ 



Jo dc ° s9 » d q 2 dcos8, + Ll dc0s9 »W 



d 2 B 



(4.5) 



d cos 9 U 



where 9^ is the angle between the /x + and the B° in the dimuon center-of-mass frame. 
We can write App m the form 



A FB {q 2 



N(z) 
dB/dz 



(4.6) 



where the numerator is given by 



N(z) 



B 



Nsm + N sm-VA + N V a + Nsp-t 



(4.7) 



The terms suppressed by m^/mb and m s /rrib have been neglected as before. Again 
for the detailed expressions, we refer the reader to Appendix [A[ 

Fig. |2| shows App{q 2 ) and the DBR for B® — > X s ^i + fi~ in the presence of NP 
in the form of Rv,a couplings, which are the ones that can most influence these 
observables. Enhancement or suppression of the DBR by a factor of 2 is possible. 
The NP couplings can enhance App up to 30% at low q 2 , make it have either sign, and 
even make the zero crossing disappear altogether. At high q 2 , however, App can only 
be suppressed. The R' VA couplings can only affect these observables mildly: a 50% 
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Figure 2: The left (right) panels of the figure show Afb and DBR for B® — > X s fi + fi~ in 
the low-q 2 (high-g 2 ) region, in the scenario where only (Rv,Ra) terms are present. The 
band corresponds to the SM prediction and its uncertainties; the lines show predictions for 
some representative values of NP parameters (Rv,Ra)- For example, the blue curves in 
the low-g 2 and high-g 2 regions correspond to (—6.85,8.64) and (—9.34,8.85), respectively. 



enhancement in DBR is possible (no suppression), but Aps can only be marginally 
enhanced and a positive zero crossing in the q 2 = 2-4 GeV 2 region is maintained. The 
mild effect of B! VA couplings as compared to the Ry,A couplings is a generic feature 
for almost all observables. This may be attributed to the bounds on the magnitudes 
of these couplings: from Sec. |2T2] , while \Rv,a\ < 10> ^ ne values of \R' V a\ < 5. 

Eq. (|4.7| ) shows that if SP or T couplings are individually present, their con- 
tribution to Apb is either absent or suppressed by m^/m^. In such a case, though 
they can enhance the DBR (marginally for SP, by up to a factor of 2 for T), Afb 
is suppressed in general (marginally for SP, significantly for T). However if both SP 
and T operators are present, their interference term is not suppressed and some en- 
hancement of Apb is possible. This still is not significant, since the magnitude of 
the SP couplings is highly constrained from — > fi + measurements. A positive 
zero crossing in the low-g 2 region is always maintained. This may be seen in Fig. |3|. 



4.2 Polarization fractions f L and f T 

In Ref. |$3J it was pointed out that, besides the dilepton invariant mass spectrum and 
the forward-backward asymmetry, a third observable can be obtained from B® — > 
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Figure 3: The left (right) panels of the figure show Afb and DBR for — > X s /.i + fi~ 
in the low-g 2 (high-g 2 ) region, in the scenario where both SP and T terms are present. 
The band corresponds to the SM prediction and its uncertainties; the lines show pre- 
dictions for some representative values of NP parameters (Rs, Rp, R' s , R' p ,Ct,Cte)- 
For example, the magenta curves in the low-q 2 and high-q 2 regions correspond to 
(-1.23, -1.79, -0.86, -1.85, 0.27, -0.36) and (-1.23, -0.23, -1.35, 0.08, 1.37, 0.01), respec- 
tively. 



X s fi + fi , namely the double differential decay width: 
d 2 B 3 



dz d cos 9/j, 



[l + cos^^HTiz) + 2 cos 6 ^HAiz) + 2(1 -cos 2 6^) H L (z) . (4. 



The functions Hi(z) do not depend on cos 6^. The sum Hl(z) + Ht(z) gives the 
differential branching ratio dB/dz, while the forward-backward asymmetry is given 
by 3Ha/4:(H l + H T ). Splitting dB/dz into longitudinal and transverse parts sepa- 
rates the contributions with different q 2 dependences, providing a third independent 
observable. This does not require measuring any additional kinematical variable - 
q 2 and cos6* M are sufficient. Including all the NP interactions, and neglecting terms 
suppressed by m^/mt, and m s /mb, Hl(z) and Ht(z) are given by 

H L (z) = H s L M (z) + H[ M - VA (z) + Hl\z) + H s L p (z) + H T L (z) , (4.9) 

H T (z) = H* M (z) + H* M - yA (z) + H%\z) + H* p (z) + H%(z) , (4.10) 

where the H functions are given in Appendix |A]. The superscripts indicate the 
Lorentz structures contributing to the term. The polarization fractions fi and fx 
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Figure 4: The left (right) panels of the figure show fi for — > X s fj, + fi~ in the low- 
q 2 (high-g 2 ) region, in the scenario where only (Rv,Ra) terms are present. The band 
corresponds to the SM prediction and its uncertainties; the lines show predictions for some 
representative values of NP parameters (Rv,Ra)- For example, the blue curves in the 
low-q 2 and high-g 2 regions correspond to (—8.14,5.75) and (1.87,4.85), respectively. 



can be defined as 

/l = _^W_ }T = glif) (41!) 

U H L (z) + H T (z) ' ,T H L (z) + H T (z) 1 > 

In the SM, fi can be as large as 0.9 at low q 2 , and it decreases to about 0.3 at high 

Fig. |] shows that when only Ry,A couplings are present, in the low-g 2 region fi 
can be suppressed substantially, or even enhanced up to 1. A similar effect - small 
enhancement or a factor of two suppression - is possible at high q 2 . The suppression 
at low-g 2 is typically correlated with an enhancement at high-g 2 . The effect of R' v A 
couplings is similar, but much milder, as expected. SP and T operators, individually 
or together, can only have an marginal effect on f^. 



5. B° 



In this section we examine the NP contributions to the radiative leptonic decay 
B® — > ji^^f. This decay has not been detected as yet. The SM prediction for the 
BR in the range g 2 < 9.5 GeV 2 and g 2 > 15.9 GeV 2 is w 18.9 x 10~ 9 [gTJ. Although 
this decay needs the emission of an additional photon as compared to B° s — > , 
which would suppress the BR by a factor of a em , the photon emission also frees it 
from helicity suppression, making its BR much larger than — > . 



This decay has contributions from many channels [|3§ 



39 



real 



direct emission of real or virtual photons from valence quarks of the B®, (ii 
photon emitted from an internal line of the b — > s loop, (iii) weak annihilation due 
to the axial anomaly, and (iv) bremsstrahlung from leptons in the final state. The 
photon emission from the b — > s loop is suppressed by m^/m 2 ^ [[37 1, an d the weak 
annihilation is further suppressed by AqcD/ m b |^T]- These two contributions can 
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then be neglected. The bremsstrahlung contribution is suppressed by m^/m^, and 
dominates only at extremely low photon energies due to the infrared divergence. The 
virtual photon emission dominates in the low-g 2 region around the resonance. If 
we choose the regions 2 GeV 2 < q 2 < 6 GeV 2 and 14.4 GeV 2 < q 2 < 25 GeV 2 as 
the low-g 2 and high-g 2 regions, respectively, then the dominating contribution comes 
from the diagrams in which the final-state photon is emitted either from the b or 
the s quark. Then the B® — > // + /i~7 decay is governed by the effective Hamiltonian 
describing the b — > sfi + u,~ transition, as given in Eq. ( p.l|) , and our formalism is 
applicable. Here we consider the the DBR and A FB in — )■ fi + fi~j. 



5.1 Differential branching ratio and forward-backward asymmetry 

We begin with the differential branching ratio. The SP operators do not contribute 
to the amplitude of B® — > /U + //~7 and hence do not play any role in the decay. 

In terms of the dimensionless parameter x 7 = 2E 1 /rriB s , where E 1 is the photon 
energy in the B° s rest frame, one can calculate the double differential decay rate to 
be 



d 2 r 



dx 7 d(cos6' /J ) 2rriB s (8tt 



1 2vm k^ M: 



(5.1) 



where v 



1 — 4m 2 /[m 2 3s (l — x 7 )]. From Eq. ( |5.1| ) we get the DBR to be 
dB 

dX-y 



1 ^ dcos# 
1 dx 7 d(cos6' At ) ' 



T B 3 



2vm 2 B 



2m Ba (8tt) 



1 1QGI a 2 , 



2 16tt 2 



em |t/ i/*|2 2 
VtbVtsl e 



Here the quantity has the form 



. 



(5.2) 



= - m 



XvA + Xt + XvA-T 



(5.3) 



where the X terms are given in Appendix [B]. The subscripts of the X terms denote 
the Lorentz structure (s) contributing to that term. For the sake of brevity, we have 
included the SM contributions in Xva- 

The normalized forward-backward asymmetry of muons in B° s — > u J + u J ~ / y is de- 
fined as 



where 9^ is the angle between the three-momentum vectors of the B® and the fi + in 
the dimuon center-of-mass frame. The calculation of Afb gives 



A FB {q 2 ) = q ( 2m B 



YvA + YvA-T 



(5.5) 



with the Y terms are defined in Appendix ||. 

The details of the calculation are given in Appendix [FJ. For the numerical cal- 
culations, we use the matrix elements given in Ref. fiOfl . The parameters involved 
in the form factor calculations are chosen in such a way that the LEET relations 
between form factors are satisfied to a 10% accuracy [|||. In our numerical analysis 
we take the errors in these form factors to be ±10%. 

Within the SM, Ap B (q 2 ) is predicted to vanish around q 2 ~ 4.3 GeV 2 (i.e. 
x.y ~ 0.85) [f4(J, and the crossing is predicted to be negative. It is therefore interesting 
to see the effects of various NP operators and their combinations on Afb- I n the 
extreme LEET limit, using the form-factor relations given in Ref. [ID], one can easily 
see that the Afb is independent of the form factors. In Fig. ||] we see large bands in 
the SM predictions of Afb in the low q 2 region. One may tend to interpret these as 
large corrections to the LEET limit, however this would be somewhat misleading, as 
we take the errors in the form factors, due to corrections from the LEET limit, to 
be uncorrelated. In realistic models, LEET corrections to the form factors will be 
correlated, leading to a smaller uncertainty band for Afb m the SM. 

Fig. H] also shows Apb and DBR in the presence of NP in the form of Rv,a 
couplings. With the large allowed values of \Rv,a\ and the absence of any helicity 
suppression, we expect VA operators to have a significant impact on the observables. 
As can be seen from the figure, the maximum allowed value of DBR can be 2-3 
times larger than the SM prediction. The BR can also be suppressed below the SM 
prediction due to destructive interference. In the low-g 2 region, the suppression can 
be large. The features of the zero-crossing predicted by the SM can be affected: it 
can be positive or negative, can take place at any value of q 2 , and can disappear 
altogether. As expected, the impact of R' V a couplings is much milder. In particular, 
the zero-crossing is always positive and in the low-g 2 region. 

With new tensor couplings, an enhancement of the DBR by up to a factor of 3 
in comparison to the SM prediction is possible. Moreover, in the limit of neglecting 
the muon mass, T operators do not contribute to the V-terms in Eq. ( |5.5[ ); their 
contribution is only to 0. As a result, they can only suppress Afb from its SM 
value. 

When all NP operators are allowed, we find that B(B® — > /! + /i~7) can be en- 
hanced by a factor of 4, or it can be suppressed significantly. The shape of Apsiq 2 ) 
is determined by the new VA couplings, while its magnitude can be suppressed if the 
T couplings are significant. 
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Figure 5: The left (right) panels of the figure show Aps and DBR for — > /i + /x~7 in the 
low-g 2 (high-g 2 ) region, in the scenario where only (Rv,Ra) terms are present. Note that 
here q 2 = w^{l — x 7 ). The band corresponds to the SM prediction and its uncertainties; 
the lines show predictions for some representative values of NP parameters (Rv, Ra)- F° r 
example, the magenta curves in the low-g 2 and high-q 2 regions correspond to (2.47, 7.08) 
and (—7.14,-0.42), respectively. 



6. B° d KjjL + ji~ 

The decay mode B d — > K{i + fi~ is interesting primarily because the forward-backward 
asymmetry of muons is predicted to vanish in the SM. This is due to the fact that the 
hadronic matrix element for the B® — >• K transition does not have any axial- vector 
contribution. Afb can have a nonzero value only if it receives a contribution from 
new physics in the form of SP or T operators. Thus, the information from this decay 
is complementary to that from the other decays considered earlier, which were more 
sensitive to new physics VA operators. 

The total branching ratio of — > Kfi + [i~ has been measured to be [75 



B(B° d /?/iV) = (4.511; 2 ,) x irr 7 , (6.1) 

which is consistent with the SM prediction |32 



B(B° d -> Kn + H~)su = (3.5 ± 1.2) x 10" 7 . (6.2) 



The integrated asymmetry, (Afb), has been measured by BaBar [0] and Belle [10 



85] to be 

(A FB ) = (0.15+°; 2 ^ ± 0.08) (BaBar) , (6.3) 
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(A FB ) = (0.10 ± 0.14 ± 0.01) (Belle). 



(6.4) 



These measurements are consistent with zero. However, within 2a they can be 
as large as ~ 40%. Experiments such as the LHC or a future Super-5 factory will 
increase the statistics by more than two orders of magnitude. For example, at ATLAS 
at the LHC, after analysis cuts the number of B® — > Kji + /i~ events is expected to be 
~ 4000 with 30 fb _1 of data |86|| . Thus, (Afb) can soon be probed to values as low 
as 5%. With higher statistics, one will even be able to measure Afb as a function of 
the invariant dimuon mass squared q 2 . This can provide a stronger handle on this 
quantity than just its average value (A F b)- 

The effect of NP on (Afb) and the Afb{q 2 ) distribution in B® — > Kfi + was 
studied in Refs. |^8[ and [[19j respectively. In the latter, it was shown that simulta- 
neous new-physics SP and T operators can lead to a large enhancement of Afb{q 2 ) 
in the high-g 2 region. However, NP effects due to other operators were not studied. 
Here we present a complete analysis of the effect of NP on the ApB^q 2 ) distribu- 
tion in B® — > Kfi + fi~ by taking into account all possible NP operators and their 
combinations. In addition, we study the possible zero crossing of ApB{q 2 ) and the 
correlations between the DBR and Afb features. 



6.1 Differential branching ratio and forward-backward asymmetry 

The differential branching ratio for this mode is given by 



dB 

dz 



XyA + ^SP + X' T + X[ 



VA-SP 



VA-T 



(6.5) 



where the normalization factor B' , the phase factor <p and the X' terms are given 
in Appendix |C| The subscripts for the X' terms denote the Lorentz structure(s) 
contributing to that term. 

The normalized forward-backward asymmetry for the muons in B® — > K[i + 
is defined as 



A FB {q 2 



l dcosd d*B 
d dq 2 d cos 6fj 



o 



d 2 B 



x ^ dq 2 d cos 6 



ii 



l dcose d 2 B 
u M dq 2 d cos Op 



° dcos9 d 2 B 
_i M dq 2 d cos 9^ 



(6.6) 



where 9^ is the angle between the three-momenta of the B® and the fi + in the dimuon 
center-of-mass frame. The calculation of ApB(q 2 ) gives 



A FB {q 2 



2-Bq 0ti ' 



J- M 



dB/dz 

where the Y terms are given in Appendix 



VA-SP 1 VA-T 



SP-T 



(6.7) 
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The largest source of uncertainty in the calculations are the B — >• K form factors. 
As these cannot be calculated from first principles within QCD, one has to rely on 
models. In the numerical calculations, we use the form factors as calculated in 
Ref. in the framework of QCD light-cone sum rules; the details are given in 
Appendix |CJ There are, however, certain limits in which relations between form 
factors can be rigorously obtained. In the large energy (LEET) limit, these relations 
are valid up to a s , 1/E K and l/m b corrections [E0[ ^T| . 

In the LEET limit, using the form-factor relations in Eq. (|C.10 ), one can verify 
that Aps is independent of the form factors. This is quite useful as it implies that 
the measurement of Afb can be used to extract the parameters of the new-physics 
operators without form-factor uncertainties in this limit. 

In the low-energy, large q 2 , region one can also derive relations between form 



factors in the heavy-quark limit [[F3|, |74| . However, these relations do not completely 
eliminate the form-factor dependence of the calculated quantities, and hence we do 
not consider these relations. An analysis where these relations have been used in the 



context of b — > sji + jj, can be found in Refs. p| 



From Eq. (|6.?1 ), clearly new VA couplings alone cannot give rise to Afb, which 
vanishes in the SM in any case. Note that this is one of the few cases where the VA 
couplings fail to influence an asymmetry significantly, in spite of the large allowed 
values of the couplings. This is because the argument about the hadronic matrix 
element B® — >■ K not having any axial-vector contribution stays valid even in the 
presence of NP. The DBR can, however, be enhanced by up to a factor of 2, or 
marginally suppressed. 

The contribution of SP operators through the Y VA _ SP terms can give rise to Afb, 
where the VA contribution comes from the SM operators. The effect is rather small 
when only Rs,p or only R' s P couplings are present, due to the strong constraints on 
their values. The peak value of Afb in the low-g 2 region stays below the percent 
level, while in the the high-g 2 region it can be enhanced up to 2% at the extreme 
end point (q 2 > 22 GeV 2 ), which is virtually impossible to observe. However if both 
the primed and unprimed SP couplings are present simultaneously, the constraints 
on them are weakened. In such a situation, the peak value of Afb in the low-g 2 
(high-g 2 ) can become ~ 5% (~ 3%). This may be seen in Fig. It is also observed 
that Afb is always positive or always negative, i.e. there is no zero crossing. The 
DBR also is significantly affected only if both the primed and unprimed SP couplings 
are present: it can be enhanced by up to a factor of 3. 

New T couplings are also expected to give rise to A F b through the Yy A _ T terms 
in Eq. (p.!2|) . It is observed from Fig. [F]that A F b(<1 2 ) can be enhanced up to 5-6% 
in almost the entire q 2 region. Moreover, at q 2 ^ 21 GeV 2 , the peak value of ApBiq 2 ) 
reaches a larger value ( ~ 30%). The value of ApB^q 2 ) is always positive or always 
negative, i.e. there is no zero crossing point. The DBR values do not go significantly 
outside the SM-allowed range. 
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Figure 6: The left (right) panels of the figure show App and DBR for — >• Kfi + fi~ in 
the low-g 2 (high-g 2 ) region, in the scenario where all NP SP couplings are present. The 
band corresponds to the SM prediction and its uncertainties; the lines show predictions 
for some representative values of NP parameters (Rs , Rp , R' s , R' P ) ■ For example, the 
blue curves in the low-g 2 and high-g 2 regions correspond to (—2.50, 6.18, —2.84, —5.64) and 
(-2.41, 1.86, -2.07, 1.42), respectively. 



When VA and T couplings are present simultaneously, a DBR enhancement of up 
to a factor of 2 is possible, while App can be large only at extremely high q 2 . On the 
other hand, when SP and T couplings are present simultaneously, their interference 
can have a large impact on App. The interference term Y' SP . T that contributes to 
App is not suppressed by m^/rrib, and therefore a large App is possible, as can be 
seen from Fig. This is also the only combination of NP couplings where a zero 
crossing may occur. Among the asymmetries considered in this paper, this is the 
one where the SP and T operators can have the largest impact. The DBR can also 
be enhanced by up to a factor of 2-3 at large q 2 due to the simultaneous presence of 
primed and unprimed SP operators. 



7. B* 



+ ,.- 



The measurement of the forward-backward asymmetry in — > K*fj, + fi by the Belle 
collaboration [110, 11 1, which showed a deviation from the SM prediction, indicates 



the possibility of the presence of new physics. According to the SM, App is < 20% 
and negative at low q 2 , has a zero crossing at q 2 pa 4 GeV 2 , and is positive but 
< 40% for larger q 2 values. The experiment showed the asymmetry to be positive 
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Figure 7: The left (right) panels of the figure show Afb and DBR for — > K/j, + fi~ 
in the low-q 2 (high-g 2 ) region, in the scenario where only T terms are present. The band 
corresponds to the SM prediction and its uncertainties; the lines show predictions for some 
representative values of NP parameters (Ct,Cte)- For example, the blue curves in the 
low-g 2 and high-g 2 regions correspond to (0.30,0.37) and (0.49,0.57), respectively. 



throughout the range of q 2 - consequently no zero crossing - and Afb 
large q 2 values. This has generated a special interest in this decay. 



60% at 



There have already been a number of theoretical studies, both within the SM [52 



53. 60] and in specific NP scenarios [58], ^HJ 32], focusing on the branching fraction 
and Aps of B\ —> K*n + ji~ . For example, Ref. [^9| has pointed out that Afb(q 2 



is 



a sensitive probe of NP that affects the SM Wilson coefficients. Other observables 
based on the K* spin amplitudes of this decay are at present under active theoretical 
and experimental analysis |57|, [58], |60|. Finally, more challenging observables, such 
as the polarized lepton forward-backward asymmetry [46, [|7], |55], E(J, have also been 



considered, though the measurement of this quantity is still lacking. 

In the coming years, the LHCb experiment will collect around 3000 events of 
B® — > K*[i + fi~ per fb _1 in the full range of q 2 . An integrated luminosity of 2 fb _1 
already would allow the extraction of the SM zero of Afb (if it is there) with a 
precision of ±0.5 GeV 2 j87|. Indeed, a dataset of 100 pb -1 would already improve 
the world precision obtained by Babar, Belle and CDF. These measurements would 
also permit many of the additional tests for NP mentioned above. 

The decay B® — > K*fi + {i~, with K* decaying to Kir, has four particles in the 
final state. This implies that there are three physical angles that can specify the 
relative directions of these four final-state particles. The differential decay rate as 
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Figure 8: The left (right) panels of the figure show Afb and DBR for — > K/j, + fi~ 
in the low-g 2 (high-g 2 ) region, in the scenario where both SP and T terms are present. 
The band corresponds to the SM prediction and its uncertainties; the lines show pre- 
dictions for some representative values of NP parameters (Rs, Rp, R' s , R' p ,Ct,Cte)- 
For example, the magenta curves in the low-q 2 and high-<7 2 regions correspond to 
(-0.09, -2.24, 0.16, -2.14, -0.33, -0.40) and (-0.40, 1.87, -0.59, 1.88, -0.34, 0.66), respec- 
tively. 



a function of these three angles has much more information than just the forward- 
backward asymmetry. Indeed, Afb is just one of the observables that can be derived 
from the complete angular analysis of this decay. In this section we also consider 
other CP-conserving observables. 

7.1 Angular analysis 

The complete angular distribution in B d — > K*fi + fi~ has been calculated in Refs. [38, 



89| within the SM. In this section, we calculate the angular distribution in the pres- 
ence of NP, which is a new result. The full transition amplitude for B(ps) — > 
K*{p K »,e*)(i + {p+)n~{p-) is 



iM(B° d K*fj, + fi- 



-i) 



4 G F a en 
y/2 4tt 



iKVt 



lb) 



[M Vfl L» + M^U* + M S L + M P L 5 + M T ^ + iM Bp L a/J e^] , (7.1) 

where the L's are defined in Eq. ( |2.8j ). The M's are given in Appendix [D|. 

The complete three-angle distribution for the decay B — » K*(—> Kn)fi + fi~ can 
be expressed in terms of q 2 , two polar angles 0„, 9k, and the angle between the 
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Figure 9: The description of the angles Ou,,K and <f> in the angular distribution of B — > 
K*{-^- Kir)fj, + fi~ decay. 

planes of the dimuon and decays, 0. These angles are described in Fig. || We 
choose the momentum and polarization four-vectors of the K* meson in the dimuon 
rest frame as 

p K * = (E K *,0,0, \p K *\) , 

e(0) = —(|#H,0, , e(X = ±1) = T-^=(0, l,±i,0) , (7.2) 

m K * y/2 

with 

,2 ™2 Jl 



Ek* = B n * , \PkA = J Eh -m 2 K ,. (7.3) 
The three-angle distribution can be obtained using the helicity formalism: 

= Njr X 



<ig 2 <i cos 0^<i cos 0#<: 



| cos 2 K (ij + J° cos 20 M + 7° cos 0^) + sin 2 X (if + 7 2 T cos 20 M + ij cos 0^ 
+/J sin 2 M cos 20 + if sin 2 M sin 20 j + sin 26> A - ( l{ T sin 20 M cos 
+/ 2 LT sin 26^ sin + J^ T sin M cos + sin 0^ sin j j , 



(7.4) 

where the normalization factor Np is 

q„2 r<2 |T/*T7 l 2 !^ 



Here /3 M = Jl — Am 2 /q 2 , and i s the magnitude of the K"* momentum in the 

B- meson rest frame: 

1 



\Pk*\ = 7) \/ m B + m K* + <? 4 - 2[q 2 m 2 B + m 2 K *(m 2 B + q 2 )] . (7.6) 
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The twelve angular coefficients I depend on the couplings, kinematic variables and 
form factors, and are given in Appendix [D]. In this paper we concentrate on the 
CP-conserving observables: the DBR, the forward-backward asymmetry Ap B , the 

(2) 

polarization fraction fi, and the asymmetries A T and Alt- 

The theoretical predictions for the relevant B — > K* form factors are rather 
uncertain in the region (7 GeV 2 < q 2 < 12 GeV 2 ) due to nearby charmed resonances. 
The predictions are relatively more robust in the lower and higher q 2 regions. We 
therefore concentrate on calculating the angular distribution in the low-g 2 (1 GeV 2 < 
q 2 < 6 GeV 2 ) and the high-g 2 (q 2 > 14.4 GeV 2 ) regions. For numerical calculations, 
we follow Ref. |62| for the form factors: in the low-g 2 region, we use the form factors 
obtained using QCD factorization, while in the high-g 2 region, we use the form factors 
calculated in the light-cone sum-rule approach. 

7.2 Differential branching ratio and forward-backward asymmetry 

The forward-backward asymmetry for the muons is defined by 

i / 2 \ ^ COS fi f_ I d COS Op 

ApB ^ q > = ~ « ; r0 JZZTn Ipf— ■ V J > 



L dcOs9 u , A r a + f i d COS a -r~^-, 77- 

JO " dg d cos ^ J — 1 A 1 dq z d cos 0^ 

It can be obtained by integrating over the two angles 9k and in Eq. ( |7.4| ). We 
obtain the double differential decay rate as 



- (l° + 7° cos 29^ + 4° cos 9^j + (if + J 2 T cos 29, 



d 2 T _ 8nN F 
dq 2 d cos 9^ 3 

+/ 3 T cos^)] . (7. 



Further integration over the angle 9^ gives the differential decay rate. The contri- 
bution of the NP operators to the differential branching ratio and forward-backward 
asymmetry of B% — > K*fi + fi~ was examined in detail in Ref. We do not repro- 
duce the analysis here, but only give the results below. 

If only Ry,A couplings are present, A FB can be enhanced at low q 2 , while keeping 
it positive, so that there is no zero crossing as indicated by the recent data jEDl |TT], 
i~3|, |i~4|j. However, an enhancement at high q 2 , also indicated by the same data, is not 



possible. On the other hand, if only R' VA couplings are present, Afb can become 
large and positive at high q 2 , but then it has to be large and negative at low q 2 . 
These couplings are therefore unable to explain the positive values of A FB at low q 2 . 
Thus, in order to reproduce the current — > K* \i~ experimental data, one needs 
both unprimed and primed NP VA operators. The NP coupling values that come 
closest to the data typically correspond to suppressed DBR at low q 2 . (See Fig. |10| .) 
But it is also possible to have a large Ap B (up to 60%) in the entire q 2 region while 
being consistent with the SM prediction for the DBR. At present, the errors on the 
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Figure 10: The left (right) panels of the figure show Apb and DBR for 5° -»• K* H + fi in 
the low-q 2 (high-g 2 ) region, in the scenario where both (Rv,Ra) an d (R'y,R' A ) terms 
are present. The band corresponds to the SM prediction and its uncertainties; the 
lines show predictions for some representative values of NP parameters (Ry, Ra, R'viR'a)- 
For example, the red curves for Afb in the low and high q 2 regions correspond to 
(—1.55,1.75,6.16,1.73) and (—5.79,1.10,0.47,-3.34), respectively. The pink curves for 
DBR in the low-g 2 and high-g 2 regions correspond to (1.96, —4.09, 4.61, 0.13). For compar- 
ison, the experimental data are also displayed in blue cross lines. 



measurements are quite large. However, if future experiments reproduce the current 
central values with greater precision, this will put important constraints on any NP 
proposed to explain the data. 

New SP couplings by themselves cannot significantly affect either the DBR or 
the Afb predictions of the SM. New T couplings in general tend to enhance DBR 
significantly, by up to a factor of 2, while not contributing any additional terms to 
the asymmetry. As a result, the magnitude of A FB is suppressed. The zero crossing 
can be anywhere in the entire q 2 range, or it may disappear altogether. However, 
whenever it is present, it is always a SM-like (positive) crossing. When SP and T 
couplings are present simultaneously, additional contributions to A FB that are not 
suppressed by m^/m B are possible. As a result, A FB obtained with this combination 
can be marginally enhanced as compared to the case with only T operators. It is 
then possible to have no zero crossing. However, the magnitude of A BB cannot be 
large in the high-g 2 region. 
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Figure 11: The left (right) panel of the figure shows fa for — > K*fi + [i~ in the 
low-g 2 (high-g 2 ) region, in the scenario where both (Rv,Ra) and (R' V ,R' A ) terms are 
present. For example, the blue curves in the low-g 2 and high-g 2 regions correspond to 
(1.64,-0.90,4.27,-0.91) and (1.96,-4.09,4.61,0.13), respectively. For comparison, the 
experimental data are also displayed in blue cross lines. 



7.3 Polarization fraction f L 

The differential decay rate and K* polarization fractions can be found by integrating 
over the three angles in Eq. (|7.4|) to get 

dT 8nN F .... n . 

^ = ^— (A L + A T ), (7.9) 

where the longitudinal and transverse polarization amplitudes Ai and At are ob- 
tained from Eq. ( [7.8|) : 

A L ={n-\n), A T = 2(ll-\lT). (7.10) 



It can be seen from the expressions for the J's in Appendix [D] [see Eq. ( D.15| )1 that SP 



couplings cannot affect At- The longitudinal and transverse polarization fractions, 
ft and fx, respectively, are defined as 

fL = A^TT T ' fT= A^TT T - (7 ' n) 

In the SM, fi can be as large as 0.9 at low q 2 , and it decreases to about 0.3 at high 
q 2 . As can be seen from Fig. ITTl new VA couplings can suppress fi substantially: it 



can almost vanish in some allowed parameter range. 

New SP couplings cannot change the value of fi outside the range allowed by 
the SM. This may be attributed to the strong constraints on the values of these 
couplings. New T couplings tend to suppress except at q 2 ~ 1-2 GeV 2 , where the 



value of fi cannot be less than 0.5 as may be seen from Fig. 12. Since both VA and 



T couplings tend to suppress fi, their combined effect results in a similar behavior. 
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Figure 12: The left (right) panel of the figure shows fi for B® — > K* [± + [i~ in the low- 
q 2 (high-g 2 ) region, in the scenario where only new T couplings are present. The band 
corresponds to the SM prediction and its uncertainties; the lines show predictions for some 
representative values of NP parameters (Ct,Cte)- For example, the red curves in the 
low-g 2 and high-g 2 regions correspond to (0.66,-0.14) and (0.3,-0.46), respectively. 

7.4 Angular asymmetries A T and Alt 

(2) 

In this subsection we consider the two angular asymmetries A T and Alt- The first 



quantity was discussed before in Ref. p7| , while Alt is introduced here for the first 
time. 

The CP-conserving transverse asymmetry A T can be defined through the double 
differential decay rate 

d 2 T 1 dT 



dq 2 d(j) 2tt dq 2 



1 + f T (a^ cos 2(f) + A^ sin 2^) . (7.12) 



Here A^ 171 ^ depends on the imaginary part of a certain combination of amplitudes 
and can be used to construct CP-violating observables. We will not consider it any 
further in this work. The asymmetry A T can be obtained by integrating over the 
two polar angles 6^ and 9k in Eq. ( [7.4|) . It can be expressed as 



4* = (7.13) 

i<y\ 

We observe that A T cannot be affected by SP couplings. 
In the SM, 

(2) 4/3 2 (K| 2 -|4| 2 + |^| 2 -|Afp 

A? « -V * ■ (7-14) 



The transversity amplitudes Au± are defined through Eqs. ( p.ll|) and (|D.12j ) given in 



Appendix 0. At leading order in Aqcd/Ek*, A-QCD/'mb and a s (the LEET limit), one 
can use the form-factor relations of Refs. |50| , 5jJ and neglect terms of 0(m 2 K ,/m 2 B ) 
to obtain 

w , A\ w . (7.15) 
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Thus, in the low-g region, 

A\ « ^= , A\^-^= for i = V,A, (7.16) 

(2) 

which corresponds to the LEET limit. A T w in the SM and is independent 
of form factors up to corrections of order Aq C d/Ek*, ^QCD/^b and a s , i.e. the 
hadronic uncertainty is small. This can be seen in Figs. 13 and 14. This indicates 

(2) 

that corrections to the LEET limit are small, and makes A T an excellent observable 



to look for new-physics effects [57 



We now examine the longitudinal-transverse asymmetry Alt, defined by 



, ^— 7T/z ' w u " uq~u,(pu coa ox u — -L " aq^d(J)d COS ' /— 1 — \ 

One can compare to A^b. In the angle is integrated over its entire 
range, while in Alt is only integrated over the range (— tt/2, tt/2). This choice 
of integration range eliminates all terms which depend on the imaginary part of 
combinations of amplitudes in the angular distribution. (These eliminated terms 
can be used to construct CP- violating observables and will not be discussed here.) 
In Alt only the CP-conserving parts of the angular distribution survive. Note that, 
in the CP-conserving limit, Alt is the same as the observable S5 defined in Ref. |HJ , 
apart from a normalization constant. The quantity Alt can also be expressed in 
terms of the observables AjP and A^p defined in Ref. |60|| . However, A L t is easily 
extracted from the angular distribution and has different properties in the LEET 
limit than Aip and AjT. 

Using Eq. (|7.4|), the asymmetry A L t can be expressed as 

tLT 

We observe from Eq. ( |D.16| ) that Alt depends on the VA couplings, as well as on 
V-S, S-TE, and P-T interference terms. In the SM, 

ALT V2(A L + A T ) ' {1 - 19) 

Now, in the LEET limit, Ay A w 0. Hence, in this limit, 

a LEET „ Re[A^y4 A * + A^A^*] 

lt 00 1 — —1 • l'- zu J 

Al + At 

From this it can be shown that the SM predicts Alt = at 

g 2 w g 7 e// m b m| ^ x g6 ^ (7 

Cj nib + Cg tub 
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Thus, just like Afb, the quantity Alt also has a zero crossing which is independent 
of form factors in the LEET limit. Note that the zero crossing of Alt is different 
from that of Afb- Figs. [13] and [TJ] also demonstrate that the zero crossing of Alt 
has a very small hadronic uncertainty. This indicates small corrections to the LEET 
limit, making the position of the zero crossing of Alt a robust prediction of the SM. 
This quantity would therefore be very useful in searching for new-physics effects. 

New VA couplings can affect A T significantly: they can enhance its magnitude 
by a large amount, change its sign, and change its g 2 -dependence. The zero-crossing 
point may be at a value of q 2 different from that predicted by the SM. 

Since Alt here is identical to the observable 5*5 in Refs. 0, |67j in the CP- 
conserving limit (apart from a normalization factor), the zero-crossing in both of 
these observables is expected to take place at the same q 2 . Indeed, the results agree 
at LO, while the NLO corrections can shift the q 2 at the zero-crossing to q 2 = 2.241q os 
|3T|j . Note that the deviation due to new VA couplings can be much larger than the 



effects due to NLO corrections. 

Except at very low q 2 , the magnitude of A L t is generally suppressed by new 
VA couplings. The primed VA couplings can be constrained by Alt better than the 
unprimed VA couplings. In both cases, the value of Alt can be anywhere in the q 2 
range, and can be positive or negative. In particular, there may or may not be a 
zero crossing, and if there is, its position can be different from that of the SM. 

(2) 

New SP couplings do not affect A T , and Alt qualitatively behaves similarly 
to the SM. New T couplings in general tend to suppress the magnitudes of both 
asymmetries (see Fig. |H|) . 



8. Discussion and Summary 

Flavor-changing neutral current (FCNC) processes are expected to be incisive probes 
of new physics. In the SM, they occur only at loop level, and hence are suppressed. 
This may allow the new-physics (NP) effects to be identifiable. Of course, since 
we have no clue about what form the NP takes, the observations from a variety of 
processes are necessary. In this paper, we have focussed on the processes that involve 
the effective transition b — > sli + . 

The transition b — > is responsible for many decay modes such as B® — > 

/J> + Lt~, B Q d -> X s /j, + h~, B° s -> /i + /i~7, B Q d -> A> + /i", B Q d -»■ K*/j, + /j,~. While some of 
these processes (e.g. B° s — > fi + lT) have not yet been observed, the upper bounds on 
their branching ratios have already yielded strong constraints on NP. Some of these 
processes have been observed and the measurements of their branching fractions, 
as well as of additional observables such as the forward-backward asymmetries, are 
available. Indeed, the recently-observed muon forward-backward asymmetry in B d — > 
K*fi + fi~ has been found to deviate slightly from the SM predictions. If this is in 
fact due to the presence of NP, such NP should contribute to all the other decays 
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Only K v ,a and R v A present Only R v A and Ry A present 




« ! (Gev ! ) 9 2 (Gev 2 ) 

Figure 13: The left (right) panels of the figure show A^ and Alt for — > K* n + ii~ in 
the \ow-q 2 (high-g 2 ) region, in the scenario where both (Rv,Ra) and (R'y,R' A ) terms 
are present. The band corresponds to the SM prediction and its uncertainties; the 
lines show predictions for some representative values of NP parameters (Ry , Ra, R'y, R'a)- 
For example, the pink curves for A^p in the \ow-q 2 and high-g 2 regions correspond 
to (1.96,-4.09,4.61,0.13) and (1.64,-0.90,4.27,-0.91), respectively. The red curves 
for Alt in the low-q 2 and high-g 2 regions correspond to (—1.55,1.75,6.16,1.73) and 
(-5.79, 1.10, 0.47, -3.33), respectively. 



involving the effective transition b — > s[i + fi~ . The effects of this NP on these decay 
modes would be correlated, and hence a combined analysis of all these decay modes 
would be invaluable in discerning the type of NP present. 

While specific models of NP may be used and their effect on the relevant ob- 
servables studied, we have chosen to explore the NP in a model-independent way, in 
terms of the Lorentz structures of the NP operators that contribute to the effective 
b — > sfi + fi~ Hamiltonian. We have performed a general analysis that includes NP 
vector-axial vector (VA), scalar-pseudoscalar (SP), and/or tensor (T) operators. We 
have computed the effects of such NP operators, individually and in all combinations, 
on these decays. We have taken the couplings to be real and have considered the 
CP-conserving observables in this paper; the CP-violating observables are discussed 



in Ref. |72| . The aim is to find NP signals, and using them, to identify the Lorentz 
structure of the NP. As the first step towards this goal, we calculate the constraints 
on the NP couplings, and, keeping the couplings within these bounds, we look for the 
observables where the NP signal can potentially stand out above the SM background. 
It is crucial to understand this SM background, which makes it imperative to 
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« ! (Gev 2 ) 9 2 (Gev 2 ) 

Figure 14: The left (right) panels of the figure show and Alt f° r B^ K* [± + in 
the low-g 2 (high-g 2 ) region, in the scenario where only new T couplings are present. The 
band corresponds to the SM prediction and its uncertainties; the lines show predictions 
for some representative values of NP parameters (Ct, Cte)- For example, the blue curves 
for Ay in the \ow-q 2 and high-</ 2 regions correspond to (0.3,-0.46) and (—0.005,0.014), 
respectively. The red curves for Alt in the low-g 2 and high-g 2 regions correspond to 
(0.3,-0.46) and (0.66,-0.14), respectively. 

use observables whose values are predicted reasonably accurately within the SM. 
The main source of the SM uncertainties is the hadronic matrix elements, whose 
theoretical calculations often have errors of the order of tens of percent. We have 
handled this on many levels. First, we have tried to identify observables that will not 
be very sensitive to the hadronic uncertainties. For example in B® — > K[i + fi~, the 
SM prediction for the forward-backward asymmetry is simply zero, independent of 
any hadronic elements. Also, while the differential branching ratios may be strongly 
dependent on the hadronic matrix elements, the forward-backward asymmetries are 
less so. Furthermore, the large-energy effective theory (LEET) limits can be used to 
control the uncertainties in the low-g 2 region for observables like Afb and Aj! . For 
example, certain observables, such as the zero-crossing of A FB in B® — > K*jj, + jj,^ , 
can be shown to be robust under form-factor uncertainties in the LEET limit. The 
longitudinal-transverse asymmetry Alt in B® — > K*fi + fi~ also has a zero crossing in 
the SM with small hadronic uncertainties. These measurements can even be used to 
extract the parameters of the NP operators, to a very good approximation. 

Also, we focus only on the situations where the NP contribution can be so sig- 
nificant that it will stand out even if the SM errors were magnified. Our figures 
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show bands for SM predictions that include the form-factor uncertainties as quoted 
in the form-factor calculations, and these are overlaid with some examples of the 
allowed values of these observables when NP contributions are included. This allows 
the scaling of these uncertainties to be easily visualized. We identify and emphasize 
only those situations where the results with the NP can be significantly different 
from those without the NP, even if the hadronic uncertainties were actually much 
larger. Note that further inclusion of the NLO QCD corrections would affect the 
central values of the SM predictions to a small extent, while also decreasing the 
renormalization scale uncertainty. However, since our primary interest is looking for 
observables where the NP effects are large, a LO analysis is sufficient. 

Our results are summarized in Table p], for the cases where the NP has only 
one type of Lorentz structure: VA, SP or T. We note certain generic features of the 
influence of different NP Lorentz structures. 

New VA operators are the ones that influence the observables strongly in most 
cases. They typically can interfere with the SM terms constructively or destructively, 
thus enhancing or suppressing the differential branching ratios by up to factors of 2 
or 3. They also are able to enhance almost all the asymmetries, the notable exception 
being Apg in B% — > Kfi + fi~, where the VA operators cannot contribute. But for 
most other observables, this kind of NP can potentially be observed. This can be 
traced to the large magnitudes of the NP couplings still allowed by data, which in 
turn can be traced to the possibility of interference between the new VA operators 
with the SM operators that allows more freedom for the new VA couplings. Typically, 
the Rv,a couplings are constrained more weakly than the Bl v A couplings, since the 
corresponding operators have the same structure as those of the SM, allowing strong 
destructive interferences. Consequently, the operators with Ry,A couplings are more 
likely to show themselves over and above the SM background. We point out that 
the exception to this rule is the Afb in B% — > K* ji + at large q 2 , where the R' v a 
couplings can cause a larger enhancement. 

The SP operators, on the other hand, are handicapped by the stringent con- 
straints from the upper bound on B(B® — >■ fi + fi~). If only Rs,p or R' s P couplings 
are present, the constraints become even more severe. It is for this reason that, even 
when the SP contributions are unsuppressed by m^/rrib, they are not often large 
enough to stand apart from the SM background. 

The couplings of the T operators, viz. CV and Cte, are not as suppressed as 
those of the SP operators. Therefore, they typically contribute significantly to the 
DBRs. However, the interference terms of these operators with the SM operators 
often suffer from the m^/m h helicity suppression, and hence they tend to suppress 
the magnitudes of the asymmetries. 

The combination of multiple Lorentz structures in general gives rise to the com- 
bination of features of the individual Lorentz structures involved. In particular, if 
the VA operators appear in conjunction with another Lorentz structure, the effects 
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Observable 


SM 


Only new VA 


Only new SP 


Only new T 


B° s 




BR 


(3.35 ± 0.32) x 
10" 9 


• Marginal E 

• Significant S 


• Large E 

• Maximal S 


No effect 


B° d 


DBR 




• E (x2) 

• S (-=-2) 


• Marginal E 


• E (x2) 






A FB 


ZC« 3.5 GcV 2 


• E(30%) low q 2 

• ZC shift / 
disappearence 


• Marginal S 


• Marginal S 






h 


• 0.9 -> 0.3 

(low->high <? 2 ) 


• Large S at low 


• Marginal S 


• Marginal E 


B° 


-> /! V 


~7 

DBR 




• E (x2 - x3) 

• S (low q 2 ) 


No effect 


• E (x3) 






A FB 


ZC« 4.3 GeV 2 


• ZC shift / 
disappearence 


No effect 


• Large S 


d 


— > 


DBR 




• E (x2) 

• Marginal S 


• E at high q 2 


• Small effect 








Vanishes 


• No effect 


• E at low q 2 

• No ZC 


• E at high q 2 

• No ZC 


B° d 


-> iV> 


V" 
DBR 




• E (x2) 

• S (H-2) 


No effect 


• E (x2) 






A FB 


ZC« 3.9 GcV 2 


• E at low q 2 

• ZC shift / 
disappearence 


No effect 


• Significant S 

• ZC shift 






h 


• 0.9 -> 0.3 

(low— ► high tp) 


• Large S 


No effect 


• Significant S 






A ( 2 ) 


• t with q 2 

• No ZC 


• E (x2) 

• ZC possible 


No effect 


• Significant S 






Alt 


• ZC at low g 2 

• more -ve 
at large q 2 


• Significant S 

• ZC shift / 
disappearence 


No effect 


• Significant S 



Table 1: The effect of NP couplings on observables. E(xn): enhancement by up to a 
factor of n, S(-=-n): suppression by up to a factor of n, ZC: zero crossing. 
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of the VA operators typically dominate. The T operators can interfere with the SP 
operators without the m^/rrib helicity suppression, but the strong constraints on the 
SP operators hold them back. A remarkable exception is the combination of SP and 
T operators in the forward-backward asymmetry in B® — > Kfi + fi~ . This asymmetry, 
which vanishes in the SM, can be enhanced to ~ 5% at low q 2 with only SP operators, 
and can be enhanced to ~ 30% with T operators but only at q 2 ~ m 2 B . However, the 
presence of both SP and T operators allows the asymmetry to be ~ 40% in the whole 
high-g 2 region. A similar feature, though to a less-spectacular extent, is observed in 
A FB of B° d K*fi+^ |2 . 



With the large amount of data expected from the LHC experiments and B- 
factories in the coming years, we may be able to detect confirmed NP signals in 
the above processes. In that combined analysis of all these decay modes, as 

carried out in this paper, would enable us to identify the Lorentz structure of the 
NP operators. This will be important in establishing precisely what type of NP is 
present. 
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Notes added: After this paper was submitted, the CDF Collaboration reported 
the measurement of 

B(B° S -> //V) = (1.8^) x 10- 8 . (8.1) 



On the other hand, the recent LHCb update does not confirm this result | 9~If . They 
improve the present upper bound on B(B® — > to 

B(B° S -> < 1-3 x 10~ 8 (90% C.L.) (8.2) 



In addition, LHCb has measured various observables in B® — > K*/j, + [i ||92|| . Their 



measurement of the A FB distribution is consistent with the SM prediction, except 
in the high-g 2 region, where we now see a slight suppression. This is contrary to the 
measurement of Belle. That is, LHCb does not confirm the Belle result of a large 
FB asymmetry in the low-g 2 region. Thus, the jury is still out on whether NP has 
already been seen in these measurements. 
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A. Details of the B® — >• X s fj, + ^L analysis 

The differential branching ratio for B^ — > X s fi + fi~ in SM can be written as 



d7 =V 1 -*>V 1 -T X 



SM 



(2z + 1) ( ^ + l] |Cf | 2 + ( 2(1 ^ + (2* + 1) ) |C7oT 



+ (^ + 1)1^ + 12 



^ + ^ Re (cf 



(A.l) 



Here t = m^/m^ 06 and 2 = g 2 / (m^ ole ) 2 . The normalization constant -Bq is |28 



2 // , ™P°l e N2 







12 


32n 2 f(m c ) «(m c ) 




2 



(A.2) 



where m c = m p c ole /m p b ole . We use m c = 0.29 ± 0.02 [||], 5(5 X c ez/) = 0.1061 ± 



0.0017 |77j and 1^6^*1/1^1 = 0.967±0.009 g|. Here /(m c ) is the lowest-order (i.e. 
parton-model) phase-space factor in B{B — >■ X c ei/): 



f(m c ) = 1 — 8m 2 + 8m^ — — 24m^lnm c 



(A.3) 



and the function K,(rh c ) includes both the 0(a s ) QCD corrections and the leading- 
order (1/m 2 ) power correction to B{B — > X c eu) : 



/a \ 1 2a s (m 6 ) A /i(m c 
K(m r ) = 1 : g{ r mc) + 



3tt 



2m 2 



(A.4) 



Here the two functions are 



g{m c ) = (7r 2 -^)(l-m c ) 2 + | 
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h(m c ) = X 1 + — -f- [-9 + 24m 2 - 72m* + 72m<! - Ihrrl - 72mt lnm c l . (A.5) 
f{m c ) 

After including all the NP interactions, and neglecting terms suppressed by m^/m^ 
and m s /rrib, the total differential branching ratio dB/dz can be written in the form 



dB 

dF 



Total 



SM 



5. 



SM-Vi + -SyA + £>£p + B<i 



(A.6) 
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where 



B 



16 



SM-VA 

By a 
Bsp 
Bj> 



(1 - z) 2 (l + 2z) [Re(CfR* v ) + Re(C 10 R A )} 
+ 32 (1 - z) 2 Re{CfR* v ) , 
1(1- z)\l + 2z) \\R V \ 2 + \R A \ 2 + \R' V \ 2 + \R' A f 



->/ |2 



4(1- z) 2 z y\R s \ 2 + \R P \ 2 + \R' S \ 2 + \R' P 
^(1-^) 2 (1 + 2^)[|C T | 2 + 4|C T£ | 2 



(A.7) 
(A.8) 

(A.9) 
(A.10) 



Note that here we have separated the contribution of the SM VA operators (subscript 
SM-VA) from that of the NP VA operators (subscript VA), for clarity. 
The forward-backward asymmetry in B® — > X s fi + fi" is 



A FB (q 2 ) 



f Q d cos 9 



d?B 



A 1 dq 2 d cos 9fj 



d COS 9 



d 2 B 



t 1 dq 2 d cost 



Jo d cos d n d/L B os e, +I°i d cos ^ 



d?B 
d^d cos y u 



(A.ll) 



where 0^ is the angle between the fi + and the B° in the dimuon center-of-mass frame. 
We can write Aps in the form 



A FB {q 2 ) 



N(z) 
dB/dz 



(A.12) 



where the numerator is given by 



N(z) 



B 



Nsm + A SM-VA + Ny A + N S P-T 



(A.13) 



with 



Nsm 
Nsm-va 
N VA 

Nsp-T 



-8C w (l-z) 2 2Cf + zRe(Q 
-8 (1 - z) 2 zRe(c l0 R v + CfR* A J + 2CfRe{R* A ) 
-8 z (1 - ^) 2 [Re(i?yi%) - Re(i4i4*)l , 
-8^(1- ^) 2 [Re{ (R s + R P ) (C T - 2C TE )} 
+ Re\(R' s - R' p ) (C* T + 2C* TE 



(A. 14) 
(A.15) 
(A.16) 

(A.17) 



The expressions of Eqs. (|A.7|) - (|A.17|) are in agreement with Ref. f30 
The polarization fractions fp and fa are defined as 



h 



Hl{z) 



H L z)+H T (z) 



fi 



Hr(z) 



H L (z) + H T (z) 



(A.18) 
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where 



H L {z) = H[ M (z) + H[ M - VA (z) + H v L A {z) + H[ p (z) + H T L (z) , 



H T (z) = H* M (z) + H* M - VA (z) + H v T \z) + H* p (z) + H T T (z) . 



The components of Hl and Ht functions are 



H S L M (z 
H^ M (z 



H 



SM-VA i 



H 



SM-VA 



Hl\z 



H v T \z 



Hl p (z 



H* p (z 



H T L (z 



Hl(z 



SB' 



16 B' 



G -(i-z) 2 \cf + 2Cf\ 2 + \a 



10 1 



16 B' 



°z(i- z y 



Cf + -c c 7 

y z ■ 



<eff 



+ \c 



10 



3 

32 B' 



^(1 - z) 2 Re (Cf R v + C 10 R* A ) + 2Re(Cf R* v ) 



3 

8 5' 



^(1 - z) 2 zRe (Cf 5^ + C 10 i%) + 2Re(Cf 5^ 



3 

16 5' 



Hi-zf 



Ry\ + |-RaI + \R'v\ \R' 



3 

4 5' 



/ |2 



3 
8 5' 



I-R5I + I-Rpl + \R'q\ + l-R'r 



-^z{\ - zf [\R S \ 2 + \R P \ 2 + \R' S \ 2 + |i£| 2 
64 5' 



3 

128 5' 



°(2-z)(l-z) 2 [\C T \ 2 + 4\C TE \ 
(l-z) 2 [\C T \ 2 + 4\C TE \ 2 ] . 



(A.19) 
(A.20) 

(A.21) 
(A.22) 

,(A.23) 
(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 
(A.29) 
(A.30) 



B. Details of the B® — >• fi + fi 7 analysis 

The transition amplitude for 5j? — > // + //~7 is 

1 



^(5 s 4/iV 7 ) = H); 



4GF OL em /-rr*Tr \ 



X 



( 7 (fc)|s 7/1 &|5°(p B )> [(Cf + 5 y + + (C 10 + 5 A + 5 A )L" 5 ] 

+ (t(*0I Ws&I^Cpb)) [-(Cf + 5 y - # V )L" - (C 10 + 5 A - i4)L" 5 ] 



cff 



C e 

^eff 

+ (l{k)\^,ul,q v b\B% B )) [-2m b -^L»] 

+ (-f(k)\sa^b\B° s (p B )) [2C T IT + 2iC TE e^L a p] 



(B.l) 



-37- 



Parameter 


fv 


frv 


fA 


Jta 


^(GeV- 1 ) 


0.28 


0.30 


0.26 


0.33 


A(GeV) 


0.04 


0.04 


0.30 


0.30 



Table 2: The parameters for B® — > 7 form factors, as defined in Eq. ( |B.3| ). 



where the L's are defined in Eq. ( |2.8j ). 

In order to calculate the DBR, one needs the B° s — > 7 matrix elements and form 
factors. The matrix elements are given in Ref. |40|1 3 : 



j(k)\s lfl b\B°(p B )) = - 



6 ^[ivpa^ Q^k 



fv(q 2 ) 



(l(k)\s ltl75 b\B^p B )) = ze 



elk • q - ■ qk» 



f A (q 2 



(l{k)\suTp,<fb\B»(p B )) = ee^e^q^fMq 2 ) , 



{-l(k)\sia^q u b\B Q s (p B )) = ie 



elk-q-e* ■ qk„ 



['y(k)\scr flu b\BV(p B )) = -i e£/ivp(r 

- f T v{q 2 )£* p k° 



{frviq 2 ) - f TA (q 2 )} 



fTA(q 2 ) , 
{(q-k)e* p q° + (e* ■ q) q p k a 

(B.2) 



Here is the four-vector polarization of the photon and q = p B — k. For the 
Bg —7- / u + /i~7 form factors (i = V, A, TA, TV), we use the parameterization [40 

fB s m Bs 



fi(q 2 ) = Pi 



(B.3) 



'■Ai + 0.5m Ba (l-q 2 /m 2 Bs ) ' 

where the parameters and Aj are given in Table These values of parameters 
ensure that the large energy effective theory (LEET) relations between form factors 



are satisfied to a 10% accuracy [40|. In our numerical analysis we take the errors in 
these form factors to be ±10%. 

In terms of the dimensionless parameter x 7 = 2E~ i /m Bsl where E 1 is the photon 
energy in the B° s rest frame, one can calculate the double differential decay rate to 
be 

1 



d2r = 1 2v <^ M i M 

dx 7 d(cos^ At ) 2m Bs (87r) 3 



(B.4) 



where v = «/l — Am 2 l /[m 2 Bs (l — Xy)]. From Eq. (|B.4|) we get the DBR to be 



dB 

dx-v 



1 d 2 r 

l~Tf 7T\ d cos ^ 

1 da; 7 d(cos0 M ) 



1 2vm 2 B 



2m Bs (8 



it 



1 16G% at 



F u em 11/ ta*|2„2 

4 2 16n 2lVtbVtsl 



. 



(B.5) 



3 We use the convention e 



0123 
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Here the quantity has the form 

2 







g m B s x 7 



XvA + Xt + XyA-rj , 



(B.6) 



where the X terms are 

X VA = (|A| 2 + |£| 2 )m| s (3 - v 2 ) (1 - x 7 ) + (|C| 2 + |£>| 2 )2m|y(l - x 7 ) , 

X T = 4|£| 2 (3 - v 2 ) + 4|F| 2 m^y (1 - x^ f 

+ 16|Gf (3 - t; 2 ) + 16\H\ 2 m Ba (3 - 2t; 2 ) (1 - x 7 ) 2 
+8m|y (1 - x 7 )Re (E*F) + 32m 2 Bs (3 - 2v 2 )(l - x 7 )Re (G*if) , 
XyA-T = -24m M Re - 48m M Re (B*G) - 48m M m| s (l - x 7 )Re (B*H) . (B.7) 

Note that here, the V A subscript includes the SM operators. The parameters A-H 
are combinations of the Wilson coefficients, form factors and NP parameters, and 
are given by 



A 

B 

C 

D 
E 
F 
G 
H 



(Cg + Rv + R'v, 



fviq 2 



(Cg + R V — R'y) 



,J A {q 2 ) , 2m b Cf 



(Cio + Ra + R'a) 



fvjq 2 ) 
(r eS i t? dm I 

m B 3 

-2C T f TV (q 2 ) , 

frviq 2 ) ~ frAiq 2 ) 



+ 



frAiq 2 



= 2G 



2G 



2Cte frviq 2 ) > 

frviq 2 ) ~ frAiq 2 ) 



TE- 



(B.8) 



The normalized forward-backward asymmetry of muons in B® — > 7 is de- 
fined as 

d 2 B 



\-FB{q 



f 

Jo 



dcos9 



M dq 2 d cos ^ 



°dcos0 
1 M dq 2 d cos 9 ^ 



I 

Jo 



d cos 0, 



d 2 5 



dq 2 d cos 6* M 



+ 



°dcos0 * B 
1 ^ dq 2 d cos 6^ 



(B.9) 



where # M is the angle between the three-momentum vectors of the B° s and the fi + in 
the dimuon center-of-mass frame. The calculation of Afb gives 



A FB iq 2 ) = ^ [ 2m 4 B v x* 



YvA + YvA-T 



(B.10) 
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with the Y terms given by 

Y VA = (Re (A*D) + Re (B*C) \m Bs {\ - x 7 ) , 

Y YA _ T = -4m M [ 2Re (C*G) + 2m| (1 - a; 7 )Re (C*H) + Re (D*E) I . (B.ll) 



C. Details of the — >• K^l + ^l~ analysis 

The transition matrix element for B d — > K[i + pT is given by 

, .s 1 



4 Gj? a em 



X 



(K(p 2 ) |s 7m 6| B{ Pl )) [{Cf + R V + R' V )L» + (C 10 + + #4)^] 

+ (K(p 2 ) \sb\ B{ Pl )) [(R s + R' S )L + (R P + R' P )L 5 } 
+ (K(p 2 ) \sia^b\ B(p x )) [-2Cf(m b /q 2 )L^ 



+ (K(p 2 ) \sa^b\ B{ Pl )) [2Crir + 2iC TE e^L 



af3\ 



(C.l) 



where the L's are defined in Eq. (|2.8| ). 

The B d ^K matrix elements needed to calculate the decay rate and asymmetry 



in B d -)■ Ku+ur are [p| 



[K(p 2 ) \s^b\ B° d ( Pl )) = (2 Pl - q),f+(z) + ( 



1 - k 2 



[K(p 2 )\sia^b\B° d ( Pl )) 
(K( P2 )\sb\B° d ( Pl )) 
(K( P2 )\sa, u b\B° d ( Pl )) 



z 

(2 Pl - q)^q 2 - (m| - m 2 K )q^ 
m B (l — k 2 



q^[fo(z) - f+{z)} 
Mz) 



-fo(z) 



(2 Pl - q)^ - (2 Pl - q) u q^ 



m B + m K 



fr{z) 



m B + m K 



(C.2) 



2 /™2 



ni 



where k = m K /m B , rh b = m b /m B , q^ = ( Pl -p 2 )^ = ( P+ + P -)^, and z = q 2 / 
The form factors /+,o,T were calculated in the framework of QCD light-cone sum 
rules in Ref. [45|. The z dependence of these is parametrized by 

f(z) = f(0) exp( Cl z + c 2 z 2 + c 3 z 3 ) , (C.3) 

where the parameters f(0),c\, c 2 and C3 for each form factor are taken from Tables 
III, IV and V of Ref. Using these, the differential branching ratio is given by 



dB - R' /Al/2 R 



X VA + X' SP + X' T + X' VA _sp + X[ 



VA-SP yl VA-T 



(C.4) 
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where B' is the normalization factor: 



G 2 F a 2 T B T/ *, 2 s 
~ 2 12 tt 5 1 ** tol ' 



(C.5) 



the phase factor is 



1 + A; + 2 — 2(k 2 + + z) , 



(C.6) 



and the X' terms are given by 



X VA 



X SP 



X' 



VA-SP 



Y' 



Xrr> 



1 " \ + \Bf) + Ami \ B '\ 2 ( 2 + 2k " ~ z ) 

+ 4mlz\C'\ 2 + 8ml(l-k 2 )Re(B'C'*) , 
^T(\Ef + Pl\Df), 

^zm 2 B [?>\F'\ 2 + 2Pl(2\G'\ 2 -\F'\ 2 ) 

(1 - k 2 ) Re(B'E'*) + Re(C'E'*) , 

m B m B 

8 mfl m B <f)Re{A'F'*) . 



(C.7) 



Here m M = m^/m B and (3^ = yl — 4m^/z. The parameters A'-G" are combinations 
of the Wilson coefficients, form factors and NP parameters, and are given by 



A' = 2(C e 9 ff + R V + R' v ) f + (z) + 4C e 7 ff m b ^ 
B' = 2(C 10 + R A + R' A )f + (z), 



C = 2(Cio + Ra + R'a, 



z 

,2\ 



Mz) - AM 



^ = 2(iZ P + i^) mB(1 . k2) fo(z), 



m b 



F' = AC T 
G' = — ACte 



fr{z) 



m B (l + k) 

fr{z) 



m B {l + k) 



(C.8) 



The limits on the kinematical variables z and cos 9^ are 



-1 < cos^ < 1 , 4m; < z < (1 — k) . 



(C.9) 
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Note that in the large energy (LEET) limit, there are relations between form factors 
that are valid up to a s , 1/Ek and l/m& corrections |5(| 51]. These are 



/+(*) 

Hz) 

fr{z) 



C(m B , E P ) 



2 2 

m B — mp 



((m B ,E P ) 



I i —)((m B ,E P ) . 



mp 
m B J' 



(CIO) 



Thus, all form factors can be expressed in terms of a single universal soft form factor 
((m B ,Ep) in this limit. 

The normalized forward-backward asymmetry for the muons in B^ — > K{i + jjT 
is defined as 



A FB (q 2 



\ cose d 2 B 

11 dq 2 d cos 9, 



dcos8 



d 2 B 



i 



dq 2 d cos 9, 



o 



' dq 2 d cos 9^ 



a cos u,j, 



(C.ll) 



dq 2 d cos 9^ 



where 9^ is the angle between the three-momenta of the B® and the fi + in the dimuon 
center-of-mass frame. The calculation of Ap B (q 2 ) gives 



A FB (q 2 



2-Bq < 



dB/dz 



Yva-sp + Yva-t + Ysp-t 



(C.12) 



where 



Y' 

1 VA-SP 

V 

I VA-T 
1 SP-T 



m. 



m B 



Re(A'D' 



-Am^l - k 2 )Re(B'G'*) - izm^Re(C'G'*) 
~Re(D'F'*) - 2zRe(E'G'*) . 



(C.13) 



Note that only Yg P _ T term is unsuppressed by the muon mass. 
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D. Details of the B® —> K*fi + fj, angular analysis 



D.l Matrix elements 



The full transition amplitude for B{p B ) — > K*(p K *,e*)fi + (p~^)p (p ) is 



%M (B° d -> K*fi + p-) = (-i) 



4 GV o; erT 
4vr 



tb) 



x 



+ (K*{p K * 
+ (K*(p K * 



e)\s^b\B(p B )) [(Cf + #v + + (Cio + Ra + B! a )L^} 

e)\s ltll5 b\B(p B ))[-(Cf + R V - B! V )L» - (C 10 + - 
e)|si^g v (l + 7s)6|5(p B )) [-2Cf(m b /q 2 )L»] 
e)\sb\B(p B )} [(R s + R' S )L + (R P + R' P )L 5 } 
e)\s l5 b\B(p B )) [(R s - R' S )L + (R P - R' P )L 5 } 



\s(j^b\B{j> B )) [2C T IT + 2iC TEf r°*L af: 



(D.l) 



where the L's are defined in Eq. (|2.S|). Here q = p B — Pk* = P^ + P^- This can be 
written in the form 



4 G F q err , 
\/2 4vr 



■cos 



f6j 



X 



[M^Z/ + M Afl L^ + M 5 L + M P L 5 + M Tflu L^ + iM EfMU L aP e 



fiua/3'i 



(D.2) 



with 



Mi 



k A4 



M 



-A" + zS" £ ; + iC"e*.q(p B + ^*)m + iD"e*.q qti , 

-E"e^e* v p a K ^ + + iG"e*.q(p B + PAr .) M + iH"e*.qq» 

iS"e*.q, 
iP"e*.q, 

C T (iT['e^e* a (p B +p K ,f + z2£e^e*Y - iT£e„ ua pe*.qp%*q p ), 
C T E^Ti'e^e* a (p B +p K ^ + zT^/^V - iT^e^e* .qp a K ^\D.?>) 



The quantities A", 5", C", D",E", F", G" , S", P", and If (1=1,2,3) are related to 
the I? — >■ ^* form factors which are given below. The contribution to the transition 
amplitudes from the quantity D"(q 2 ) vanishes and that from H"(q 2 ) is suppressed 
because of the equation of motion of the muons. 
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D.2 Form factors 



The form factors for the decay amplitude for B% — > K*fi + fi [Eq. (|D.1|) 1 in terms of 
matrix elements of the quark operators are given by 



(K*(p K *,e)\s llx {l± l5 )b\B(p B )) 



_ . 2m K , 



± ie*(m B + m K *)M{q ) =F i(ps + PK*)f,e* • q 



A 3 (q 2 ) - A (q 2 ) 
A 2 (q 2 ) 



(m B + m K * 



*V A (7 

e^A^e p K ,q 



2V(q 2 ) 



* y (m B + m K *) ' 



(D.4) 



where 



A 3 (q 2 ) 



m B + m K * A , 2 m B -m K * 



2m 



K" 



-Mq 2 )~ 



2m 



K" 



A 2 (q 2 ) ■ 



(D.5) 



(K*(p K *,6)\sa^b\B(p B )) = ie^xA -T x {^y x {p B +PK-Y 



[m 



B 



m\«) 



+ 



T l {q 2 )-T 2 {q 2 )- 



Ti(g 2 )-T 2 (^)) e *V 
q 2 



[m 2 B - m\* 



W) K -qpW\ ■ (D.6) 



(K*(p K .,e)\si(T ltt ,q u (l±>y s )b\B(p B )) = 2e^ XtT e* v p x K .q tT T^q 2 ) 

± i\ e*„(m| - m 2 K «) - (p B +p K *)u,£* • q \ T 2 {q 



j (pB+PK*)nq 2 \ r . 
±ie ■ <l \ % ~ -7ZT2 — Z^T \ T m) 



\m B -m" K ^ 



(D.7) 



(K*(p K *,e)\s(l± l5 )b\B(p B )) = T 2^ e* ■ q A (q 2 ) . (D.8) 

m b 



Here we have neglected the strange-quark mass. The matrix elements are functions 
of 7 unknown form factors: A 0i i^(q 2 ), V(q 2 ), Ti j2i3 (g 2 ). 

The matrix elements Mv t A,s,p,T,E appearing in Eq. ( |Q.3| ) can be written in terms 
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of these 7 form factors, coupling constants and kinematic variables as 



A" = 
B" = 
C" = 
D" = 

E" = 

F" = 

G" = 

H" = 

S" = 

P" = 
T'l = 



rpll 

1 2 — 



rpll 

2 3 — 



2V(q 2 )(C e 9 ff + R v + R' V ) | 4m, 
m B + m K * q 



+ ^C^Uq 2 ) 



(m B + rriK^iq 2 )^ 1 + R v - R' v ) + ^C e 7 ff T 2 (q 2 )(m 2 B - m 2 K ,)] , 



A 2 {q 2 ) 



tti b + m K 

2 -^f(Cl" + Ry - i? v -)(-4 3 ( 9 2 ) - A (f)) - 2 -^Cf'T 3 ( q ') 



Q 

2V(q 2 )(C 10 + R A + R' A ) 
m B + m K * 

(m B + m^*)^i(<? 2 )(Cio + Ra- R'a) 

A 2 (q 2 



Tn B + m K 
2rriK 



"(Cio + Ra — R'a) 



-(C l0 + RA-R' A )(A 3 (q 2 )-A (q 2 ))\, 



-2(R s -R' s )^A (q 2 ) 

-2(R P -R' P )^A (q 2 ) 
-2T!(g 2 ), 

2(m|_m ^m( g 2 )-T 2 (0) 



Q 



4 / 2 2 q 2 T 3 (q 2 ) \l 

? (r 1 ( g )-T 2 ( g )- m ,_ < j j. 



(D.9) 



Also, we define 



T + = (g 2 + 2£„* v^)^' + g 2 ^' , T_ = 2|p*. I v^" • 



(D.10) 



D.3 Trans versity amplitudes 

We summarize the various transversity amplitudes that appear in the B® — > K*fi + fi~ 
angular distribution. The decay amplitude of B° d — > K*fi + fi~ depends on the K* 
polarization vector e(A) with helicity A (0, ±1). Hence, the decay amplitude can 
be decomposed into three components. Below we define the helicity amplitudes of 
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various operators with different Lorentz structures (V, A, S, P, T, TE) in Eq. (p.f 



K 


- m Ek - 

\rriK* 
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where the amplitude A vt is related to the time-like component of the virtual K* . In 
the transversity basis, the positive and negative helicity amplitudes are replaced by 
the transversity amplitudes as 



Aj = i(i4+ + 4r), ^ - l(A i 



i = V,A,T,TE. (D.12) 



2 y 1 1,1 "~ 2 
The left and right component of the transversity amplitudes of vector and axial- vector 



currents in 61 can be written as 
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Note that in the notation of Ref. [^TJ, we have the correspondence A, Q _g VA 



("sAzV-^O^o^ ±y The amplitudes As t p )V t remain the same, while As = —{^Jq 2 /N)A^ H . 
D.4 Angular coefficients 

The expressions for the twelve angular coefficients (J's) in the B® — > K*fi + fi~ angular 
distribution are summarized here according to K* helicity combinations A1A2. The 
longitudinal J°'s (AiA 2 = 00) are given by 
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The transverse I s (A1A2 = ++, , H — , — h) are given by 



2 + 



(|4| 2 + \Al\ 2 + \A(\ 2 + |4| 2 ) - 4(-2 + 01) (\A\\ 2 + |4| 2 + \Af 
+\Al E \ 2 ^j + ^(\A\\ 2 + \Al\ 2 - \A(\ 2 - \Ai\ 2 - 16(44* - A\ E Al E *)} 



yg 2 



- |4| 2 + + |4| 2 + |Ain - 4 MAjj \ 2 + \Ai\ z + LAj*f + \A 



m,. 



Re[A v ± Af + A\Af } + 4-^Re[Af (A i* - A{*) + A^(A \ h * + 4 A *)] 



'4 



(|4| 2 - \Al\ 2 + \Ai\ 2 - |4| 2 ) + 16(44* + A\ E Al E *\ 
l£ = 20 2 lm[Aj*Al + 4*4] . 

The mixed 7 LT 's (A1A2 = 0±, ±0) are given by 
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